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ABSTRACT 

We verified the off-axis jet model of X-ray flashes (XRFs) and examined a discovery of off-axis orphan 
gamma-ray burst (GRBs) afterglows. The XRF sample was selected on the basis of the following 
three factors: (1) a constraint on the lower peak energy of the prompt spectrum (2) redshift 

measurements, and (3) multi-color observations of an earlier (or brightening) phase. XRF020903 was 
the only sample selected basis of these criteria. A complete optical multi-color afterglow light curve of 
XRF020903 obtained from archived data and photometric results in literature showed an achromatic 
brightening around 0.7 days. An off-axis jet model with a large observing angle (0.21 rad, which 
is twice the jet opening half-angle, Ojet) can naturally describe the achromatic brightening and the 
prompt X-ray spectral properties. This result indicates the existence of off-axis orphan GRB afterglow 
light curves. Events with a larger viewing angle (>~ 29jet) could be discovered using an 8-m class 
telescope with wide field imagers such as Subaru Hyper-Suprime-Gam and the Large Synoptic Survey 
Telescope. 

Subject headings: stars flare-stars: gamma-ray burst: general-stars: supernovae 


1. INTRODUCTION 

Long gamma-ray bursts (GRBs) are believed to oc¬ 
cur when a very massive star dies in a highly energetic 
supernova forming a black hole and producing a rela¬ 
tivistic jet. Because of the release of a large amount of 
isotropic equivalent energy {Eiso ~ release 

in the short prompt gamma-ray phase (typically rang¬ 
ing from several seconds to several tens of seconds), the 
consideration of jet collimation of GRBs is necessary to 
explain the radiation mechanism from compact sources 
(e.g., massive stars and/or mergers). This necessity is 
supported by achromatic temporal breaks (also known 
as jet breaks) in the afterglow light curves of the GRBs. 
Ultra-relativistic collimation and a jet structure are re- 
quired to expla in the light curve temporal breaks (e.g., 
iSari et ^Il999f) . However, no direct observational evi¬ 
dence exists for this jet collimation. 

Off-axis orphan GRB afterglows are pr oduced as a nat - 
ural consequence of GRB jet production (jRhoads 111999(1 . 
The production of these afterglows is as follows: GRBs 
are collimated with rather narrow opening angles, and 
the afterglow that follows can be observed over a wider 
angular range. While the GRB and the early afterglow 
are collimated to within the original jet opening angle, 
the afterglow in the late phase can still be observed by 
an off-axis observer after the jet break. The Lorentz fac¬ 
tor, r, is a rapidly decreasing as function of time. This 
means that an observer at Oats cannot see the prompt 
gamma-ray emissions when Oots > djet but can detect 
an afterglow once T”^ equals 9obs- Here, 9jet is the jet 
opening half angle. As the typical emission frequency 
and flux decrease with time (while the jet opening half 
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angle 9jet increases with time), observers at larger view¬ 
ing angles will detect fainter afterglows at longer wave¬ 
lengths (e.g. optical and radio). Hence, the expected 
properties of GRB orphan afterglows are as follows: (I) 
prompt emissions in the high-energy band are absent, (2) 
their brightness is fainter than that of on-axis GRB op¬ 
tical afterglows, (3) they have the same optical color as 
on-axis afterglows, and (4) they show host galaxy proper¬ 
ties similar to those of on-axis GRBs. The afterglows are 
characterized by three-component light curves with ris¬ 
ing, peaking, and rapidly decaying phases. In this case, 
events intermediate between classical hard GRBs and off- 
axis orphan GRB afterglows should exist. 

A candid ate for the interm e diate events i s X-ray flashes 
(XRFs) as lYamazaki et al.l ()2002l 1200311 described in 
their off-axis jet model for explaining their nature. XRFs 
were recognized by the Wide Field Gameras (WFG, 2-28 
keV:lJager et al.llI997ll onboard the BeppoSAX satellite 
(iHeise et al.l 1200 lib The observed X-ray temporal and 
spectral properties of XRFs in the prompt phase do not 
show any differences relative to those of GRBs, except for 
the considerably lower energy values of the peak of the 
vFi, spectrum in the observer’s frame. A large number of 
XRF samples were provided by BeppoSAX and High En¬ 
ergy Transient Explorer 2 {HETE-2), both of which em¬ 
ployed wide field X-ray cameras-the WFG onboard 

BeppoSA X and the Wide-Field X-ray Monitor (WXM, 
2-25 keV iShirasaki et al.l 1200.'^ onboard HETE-2. The 
majority of HETE-2 samples (nine out of 16 XRFs) show 
a low energy of spectral peak energy < 20 keV 

((Sakamoto et al.l [200511 . The number of XRFs detected 
by HETE-2 was comparable and relatively larger than 
that of GRBs indicating that XRFs r epresent a large 
porti on of the entire GRB population ((Sakamoto et al.l 
1200511 . The observational properties of XRFs can be 
interpreted as being associated with the same phe¬ 
nomenon as classical hard GRBs and as being represen¬ 
tative of the extens ion of the GRB population to low 
peak-energy events ((Kinnen et al.ll200lt iSakamoto et al.l 
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I2005f) . To explain the aforementioned prompt observa¬ 
tional properties, three model s have been proposed for 
XRFs: a high redshift origin (lHeisel l2003D: the off-axis 
jet model (lYamazaki et al.ll2n02Ll200,‘itlZhang et al.l[2nn4 

iLamb et al.ll2005ll . which is equivalent to the unification 

scenario of AGN galaxies; and intrinsic properties (e.g., a 
subenergetic or inefficient fireba ll), which may also pro¬ 
duce on-axis orphan afterglows (iHuang et aljl2003l . 

The Swift satellite has also been detec ting many X-ray 
rich GRBs (XRRs) and XRF samples (jSakamoto et alJ 
1200811 . However, the XRF samples tend to be at the 
high-end of the distribution of BeppoSAX and 

HETE-2. This is because of the relatively higher en¬ 
ergy cove rage of the Burst Aler t Telescope (BAT, 15- 
150 keV; iBarthelmv et alJ 1200511 onboard Swift. Al- 
though the Monitorin g of All-sky X-ray Image (MAXI; 
iMatsuoka et al.l 1200911 attached to the International 
Spac e Station has been detecting XRFs (|Serino et alJ 
l20I4ll . there has been no appropriate follow-ups because 
of poor position determination. Hence, XRF studies have 
stagnated because of the lack of soft-X-ray monitoring in¬ 
struments and intensive multiwave length follow-up ob¬ 
servations. 

In this paper, we investigated the characteristics of 
XRFs on the basis of redshift measurements, multifre¬ 
quency afterglow monitoring, and spectral peak measure¬ 
ments The main objective was to verify the off- 

axis jet model and to provide feedback to ongoing and 
planned optical untargeted time- domain surveys by us- 
ing Subaru Hyper-Suprime-Gam (iMivazaki et al.l [201211 
and the Large Synoptic Survey Telescope (LSST), which 
have considerable potential for detecting off-axis orphan 
GRB afterglows. 


2. SAMPLES OF X-RAY FLASH 

We considered possible XRFs for our study and quickly 
realized XRF 020903 was the only event that has (I) 
either a measurement or an upper bound on the peak 
energy of the prompt spectrum (2) a measured 

redshift, and (3) multicolor afterglow observations ad¬ 
equately cover the early afterglow phase when achro¬ 
matic brightening of the afterglow might occur. This 
was one of the XRFs detected by HETE-2 and the hrst 
events for which an optical afterglow was detected and 
a spectroscopic redshi ft {z = 0.251) was determined 
(|Soderberg et al.l [20041) . The prompt emission had the 
lowest intrinsic spectral peak energy Epl^. of 3.31};® keV 
among all the XRF samples. Here, we employed the 
value estimated by the constrained Band fun ction with 
the 90% conhdent level ( Safemoto_et_aL I2004D . In accor¬ 
dance with the report of Sakamoto et al. (|2004[) , the light 
curve in the prompt phase exhibited a double-peak struc¬ 
ture and a lack of signals above 10 keV. Figure [T] shows 
histograms of the spectral peak energ y in the observe r 
frames detected by Swift/BAT (|Lien et al.ll2015[l . 

HETE-2 (|Sakamoto et al.l [20051) . t he Burst and Tran - 
sient Source Experiment (BATSE: iKaneko et al.l[2006D . 
the Fermi Gamma-ray Burst Monitor {Fermi/GBM; 
Ivon Kienlin et al.ll2014f l. and b'wzafcw Wide-Band All-Sky 
Monitor (WAM). The WAM E°///j^ distribution is pro¬ 
duced using values, which are available in GGN circulars 
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Fig. 1.— Distribution of for the Siui/t/BAT, HETE-2, 

the BATSE, the Fermi/GSM., and the Suzaku/WAM samples. 
XRF020903 showed the lowest j. among all five samples. 

up to Decemb er 20l43- The standar d analysis procedure 
of WAM dat a (lYamaoka et al.l 20091) are described in lit¬ 
erature (e.g., lOhno et al.ll2008t lUrata et al.ll2012all . For 
comparison with a large volume of data, we used 
instead of intrinsic spectral peak energies. The lowest 
intrinsic spectral peak energy determined from the rela¬ 
tion A^eafc=(l+^)-®peafc XRF020903 stands out from 
all GRB populations. The X-ray to the y-ray fluence ra¬ 
tio of 5.6 qualified thi s burst as an XRF. He re, we used 
the XRF definition bv iSakamoto et id] (j2005ll for HETE- 
2 events. 

Effective follow-up observations in the optical band 
were made u sing wider-field-o f -view (FOV ) instruments 
report ed by iSoderberg et al.l (|2004D and iBersier et al.l 
(|2006D . Because of a delay in the position alert and 
a large position error, the light curve sampling was 
sp arse around th e poss ible rebrightening epoch reported 
by IBersier et al.l (120061 1 . To describe the rebrightening 
epoch by using multicolor data, we added Subaru archive 
data, as described in §3. 

3. OBSERVATIONS AND DATA COLLECTION 

We collected data for the XRF 020903 afterglow 
by using Subaru archive data a nd photomet r ic re- 
sults obtained from literature (|Bersier et al.l 120061 : 

^ http://www.astro.isas.jaxa.jp/suzaku/HXD-WAM/WAM- 

GRB / results / gen. html 
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iSoderberg et al.l I2004D . To perform accurate optical 
photometry by removing contamination from the host 
galaxy, we also employed Panoramic Survey Telescope 
and Rapid Response System 1 (Pan-STARRSl) data. 
The individual data collections are summarized in the 
following subsections. 

3.1. Subaru Suprime-Cam 

The entire HETE-2 position error region was imaged 
by Suprime-Cam attached to the 8.2-m Subaru telescope. 
The Suprime-Cam camera consists of ten high-sensitivity 
2kx4k CCDs and covers a field of view of 34' x 27' 
(|Mivazaki et al.l [20021) . The first epoch of observation 
on 2002 September 3 (0.16 day after the burst) involved 
an i?c-band filter. Although the observation was per¬ 
formed under marginal airmass conditions (from 2.55 to 
2.75), the wide-held imager with the 8-m-class telescope 
provided the deep i?c-band images with three sets of 180 
s exposure. Subsequently, three color observations were 
made on the night of 2002 September 4 by using Rc—, 
Ic—, and z'-band hlters. During the night observation, 
two epochs of Ic— and z'— bands imaging and three 
epochs of i?c-band monitoring were also conducted to 
check the short-term variability of the afterglow. The 
Ic— and z'— band observations were made under the 
reasonable observing condition (e.g., in the airmass range 
from 1.34 to 1.97). The second epoch of i?c-band imag¬ 
ing was also performed under the reasonable observing 
condition (e.g., for an airmass of 1.32). By contrast, 
the hrst and third epochs of i?c-band observations were 
made under the marginal airmass conditions (2.65 — 2.22 
for the first epoch and 1.97 — 2.27 for the third epoch). 
These observations were performed with using appropri¬ 
ate dithering techniques to fill up the chip gaps in the 
camera. 

The Subaru-XMM Deep Survey (SXDS) team at¬ 
tempted to obtain deeper reference images with i?c-band 
filter on the night of October 9 and 10. However, the re¬ 
sults were marginal because of poor weather conditions 
and only two set of images taken with 180 s exposure 
on 2002 October 9 were available for scientific analysis. 
Table [T] shows the log of observations for the available 
images. Because of the limited data set, searching for 
a counterpart by using only the Suprime-Cam data was 
difficult. All raw images and related calibration data are 
available on the SMOKA (Subaru Mi taka Okayama Kiso 
Archive system [Baba. H et al. Il2002fl . 

3.2. CTIO 

We obtained R and I images taken by the wide field 
MOSAIC 11 camera on the Cerro Tololo Inter-American 
Observatory (CTIO) Blanco 4-m telescope from the 
NO AO archive system. These data were taken on 2004 
September 18 with 1800 s exposure (360s x 5) for the 
i?-band, and on 2004 September 9 with 1800 s expo¬ 
sure (450sX4) for the /-band, when was sufficiently late 
to estimate the host galaxy contribution on the Sub¬ 
aru images. The data were pro cessed using the pipeline 
of the NO AO archive system (IPierfederici et al] 12004 
iPierfederi^ 120061 1 . These two stacked images were also 
used to remove the host galaxy contamination for de¬ 
scribing t he late phase aftergl ow temporal evolution pre¬ 
sented bv iBersier et al.l (I2006D . We also attempted to ob¬ 
tain B, R and I images taken bv iBersier et al.l (|2006f) on 


2002 September 4 and 9 (~0.66 day and 5.7 day after the 
burst, respectively). However, these imaging data were 
unavailable because of a limitation of the tape reader 
device on the NOAO archive system. 

3.3. Pan-STARRSl z-band image 

Pan-STARRSl (PSl) z'-band images were also ob¬ 
tained to determine the host galaxy contamination on 
the Subaru z-band images taken in 2002 September 4. 
The PSl telescope has a 1.8-m diameter primary mir¬ 
ror, and it is located at the summit of Mt. Haleakala on 
Maui. The site and optics deliver a point-spread func¬ 
tion (PSF) with a full-width at half-maximum of about 1 
arcsec, over a seven square degree field of view. The PSl 
was used to conduct a 37r survey of the e ntire sky north 
of — 3 0° in (?'-, r'-, f'-, z ' - and i(/-band (IMaenier et al.l 
120131 : iSchlaflv et al.l 120121 : iTonrv et al.l l2012h . Because 
of the PSl 37r surveys strategy, the XRF020903 held 
was covered during the 3.5 years of survey starting from 
2010. The imag es were proces sed by the Image Pro¬ 
cessing Pipeline (lMagnieiil2006[) . and a deeper stacked 
z-band image was generated by the SWarp software 
([Bertin et al.ll200^ . The PSl Stt catalogs were also used 
to perform photometric calibration for the z-band images 
taken by Subaru. 

3.4. Late phase optical observations 

Late optical afterglow in the /?-band was also moni¬ 
tored by the 1.82 m Copernicus telescope at Mount Ekar 
(Asiago, Italy) on 2002 September 29 (26.5 days), and 
by the 3.5m Telescopio Nazionale Galileo (TNG) at La 
Palma with V-, R-, and /-band hlters on 2002 Octo¬ 
ber 2 (29.5 days). Additional BVRI images were ob¬ 
tained with the Danish 1.5-m telescope at the La Silla 
Observatory between 2002 October 10 (36.8 days) and 
14 (40.7 days) and then on October 26 (52.6 days). Dif¬ 
ferential images without host galaxy contamination were 
produced from late time images taken in 2004 through 
image convolution and by using subtraction methods of 
lAlard fc Luoton I (|1998f) . These photometric results for 
the optical afterglow agai nst with the s econdary stan- 
da rd stars from t he lis t of iHendt^ (|2002( 1 were reported 
bv iBersier et al.l (12006( 1. 

3.5. Radio observations 

Very Large Array (VLA) observations were performed 
at 8.5 GHz on 2002 September 27.22 UT and radio af¬ 
terglow was detected. Further monitoring observations 
were made with the VLA over 370 days at frequen¬ 
cies of 1.5, 4.9, 8.5, and 22.5 GHz (jSoderbersf et al.l 
l2004f) . The Very Long Baseline Array also observed 
the radio afterglow and determined the position accu¬ 
rately, as 02000 = 22^48'" 42L33912 ± 0L00003, d ^ono = 
-20°46'08”.945 ±0”.0005 ([Soderberg et a,l.l[200l . 

4. DATA ANALYSIS AND RESULTS 
4.1. Reduction 

The basic reduction of the Suprime -Cam data was per - 
formed using the SDFRED package (|Ouchi et al.ll2004ll . 
This entailed bias subtraction and flat-fielding in the Rc- 
, Ic- and z'-band by using a sky flat constructed from the 
median of the dithered science frames. After performing 
distortion correction for each object frame, we stacked 
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Fig. 2. — Rc/R, Ic/I, and z-band images of the XRF020903 field. The first column shows the Subaru images take n in the Rc (top), 
Ic (m iddle), and z-band (bottom). Reference images for removing the host galaxy component by using the algorithm of lAlard fc Lupton I 
II1998I1 are shown in the second column. The top-center, middle-center, and bottom-center panels show the reference images taken with the 
CTIO R-band, CTIO 7-band and PSl z-band filters. The third column shows Rc/R, Ic/I, and z-band differential images, respectively. 



Fig. 3.— 7J-band and 7-band transmission curves for Subaru 
Suprime-Cam (red solid line) and CTIO MOSAIC II (blue dashed 
line). The data were obtained from their respective instrumental 
Web-pages. 

the frames for each epoch and for each band path filter 
by considering the median. We performed astrometric 
calibration for the stacked images against the Stt catalog 
of PSl. 

For other data, we used the reduced data as previously 
summarized, except photometric calibration. The abso- 


TABLE 1 

Log of Subaru follow-up observations. 


Delay (Days) 

Filter 

Flux density (/i-Jy) 

0.1645 

Rc 

21.22 ± 1.09 

0.8609 

Rc 

34.87 ± 1.38 

1.0087 

Rc 

35.00 ± 1.39 

1.1378 

Rc 

29.57 ± 1.28 

35.78945 

Rc 

3.96 ± 0.46 

0.8991 

Ic 

49.21 ± 1.82 

1.0436 

Ic 

41.08 ± 1.67 

0.9569 

z 

55.67 ±2.06 

1.0955 

z 

54.10 ± 1.96 


lute photometric calibration for Rc a nd Ic wa s perfo rmed 
using standard stars from the list of lHende'iil (|2002D . For 
the photometric calibration of z-band images, we used 
the PSl Stt catalog and selected unsaturated stars around 
the afterglow location. 

4.2. Afterglow light curves 

As iBersier et al.l (1200611 reported, the host galaxy of 
XRF020903 is dominated and contaminated in afterglow 
measurements. Hence, to remove the galaxy contami¬ 
nation, we generated differential images using the HOT- 
PAN TS software, which emplo ys an algorithm presented 
by in lAlard fc Liinton I (jl998[ ). We also used special- 
purpose software based on the same algorithm and tuned 
to the Subaru/Suprime-Cam data (|Urata et al.ll2012bll 
to check the consistency. For Rc- and Jc-band data of 
Subaru/Suprime-Cam, we used the R- and /-band im- 
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Fig. 4. — Rc/R, Icjl, 2 :'-band and VLA 8.5GHz light curves of XRF020903. The dotted lines show the off-axis jet model functions 
described in §5.1. The dash-dotted line indicates the SN component reported bv IBersier et al.1112009 ). 


ages taken by the GTIO in 2004. The depth of these 
images was comparable to that of Suprime-Gam images. 
For the z-band, the reference image generated using the 
PSl Stt data was shallower than that of the Subaru im¬ 
age. However, there were the significant signals at the 
GRB afterglow position in the generated image (Figure 
El). Hence, we also generated the differential images by 
using the same code. As shown in FigureEl the quality of 
the differential images is appropriate for estimating the 
brightness of an afterglow component by the standard 
aperture photometry. There are appreciable differences 
between the transmission curves of the Suprime-Gam Rc 
and GTIO R bands. Figure |3] shows the transmission 
curves for the Suprime-Gam {Rc and Jc) and GTIO (i? 
and /) filters. This transmission difference causes a sys¬ 
tematic difference of 17% in the magnitudes for an object 
with a power-law spectrum of index —1. 

Figure m shows the optical light curves f or our Subaru 
data a nd photometric results reported by iBersier et al.l 
(!2006f): the ra d io af terglow light curve at 8.5GHz 
iSoderberg et al.l (|2004l) is also shown. There is a sig¬ 
nificant signal at 0.165 days after the burst. The use of 
only one epoch and single-color observations make it diffi¬ 
cult to describe the afterglow temporal evolution. When 
we consider a simple direct connection to the first epoch 
of R band photometry (0.660 days) made by GTIO, the 
Rc/R band light curve is flat. The Rc/R- and Jc//- band 
light curves show a clear rapid rebrightening between 0.7 
and 0.9 days. The equivalent rising power law index a 
{Fu oc t“) are ^ 2.6 in the Rc/R and ^ 2.3 in the Ic/I 
band. The multiepochs observations of the Subaru on 


2002 September 4 show a flatter and decaying light curve 
for all three bands. The Rc/R band light curve shows 
continuous decay up to ^7 days. This indicates that the 
rebrighte ning peak wa s arou nd 0.8^0.9 days after the 
burst. As iBersier et al.l (|2006[ ) demonstrated, there is an¬ 
other late-phase rebrightening peaking at 20 days in 
the R band light curve, which was interpreted as the as¬ 
sociated supernova component. The late phase i?c-band 
image taken by Subaru on 2002 October 9 also detected 
this supernova component in the differential image, and 
this component is consistent with that observed in the 
J?-band within the systematic error. 

4.3. Spectral flux distributions 

The multiband observations of XRF 020903 were used 
to determine the spectral flux distribution (SED) at 0.7 
and 0.9 day after the burst. To remove the effects of the 
Galacti c interstellar extinction , we u sed the reddening 
map of iSchlaflv fc Finkbeined (|2011l) . Because not all 
multiband observations were performed exactly at the 
same epoch, their corresponding fluxes were rescaled to 
assume a power-law function (/% cx; t~°‘). In particular, 
the brightening phase around 0.7 day was substantially 
affected. To determine the spectral energy distribution 
(SED) at 0.7 day, we fixed the time at the epoch of the 
R-band observation (0.695 days). Because the Rc/R- 
and Jc/J-band observations were made earlier than B 
band observations, and the brightening index for Rc/R 
and Ic/I were estimated as previously described (§4.2). 
In figure [5l we plot the SED obtained on the basis of 
the photometric results for BRI bands provided by the 
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Fig. 5.— Spectral flux distribution of the XRF020903 afterglow 
at 0.695 (red circle point) and 0.899 (blue square points) day. The 
dashed lines show the best-fit model functions. 

CTIO observations. The SED is well fitted by the power- 
law function as f{iy) oc , where /(i^) is the flux den¬ 
sity at frequency v, and /3 is the spectral index. We have 
obtained the /3 at 0.695 day as 1.48 ± 0.06. Similarly, 
we generated the SED at 0.899 day on basis of the on 
Rc-,Ic-, and z'-band results and obtained a j3 value of 
1.43 ± 0.08, which is consistent with the value for 0.695 
day. These results imply that the rebrightening is achro¬ 
matic. 

5. DISCUSSION 

On the basis of the results for XRF020903 described in 
the preceding section, we discuss the off-axis jet model 
in the following sections because this model provides a 
reasonable explanation for this particularly soft XRF. 
Furthermore this model can explain both the extremely 
soft prompt emission features and achromatic light curve 
brightening altogether. 

5.1. Off-axis modeling of afterglow light curves 

To perform light curves and SEDs modeling 
for XRF020903, we employed the boxfit code 
(|van Eerten et all I2012f) that involves two-dimensional 
relativistic hydrodynamical jet simulations for deter¬ 
mining the burst explosion parameters, including the 
off-axis viewing angle and the synchrotron radiation 
parameters for a homogeneous circumburst medium. 

In the modeling, we added 17% of the systematic errors 
to the i?-band light curve because of the large transmis¬ 
sion differences, as shown in Figure [3] To perform light 
curve modeling, we focused on the data between 0.6 and 
8 day because the observed optical light curve showed 
achromatic brightening between 0.7 and 0.8 day after 
the burst, which is unusual among the well-observed af¬ 
terglows. The later phase (from ^20 days) was excluded 
in the modeling bec a use it is likely to be the SN com¬ 
ponent iBersier et al.l (|2006ll . We also excluded the single 
data point at 0.16 day because the sparse monitoring and 
single detection at 0.16 day make it difficult to decode 
the light curve between 0.16 and 0.7 day. The considera¬ 
tion of different components (one component faded away 
before 0.6 day and the other was rebrightening) was rea¬ 
sonable because connecting the light curve between 0.16 
and 0.7 day smoothly produced a flat evolution, indi¬ 
cating long lasting (^0.6 days) energy injection. Using 
multicolor optical data between 0.6 day and 8 days and 



Fig. 6. — XRF080330 multi-band light curves with the off-axis jet 
mode function (dot ted lines). The optica l and near-infrared data 
points reported by IGuidorzi et al.l II2009I) were used for the light 
curve modeling. 

the radio 8.5 GHz radio data, we performed light curve 
modeling. We examined the light curves for various val¬ 
ues of 9jet and 9obs- The range of 9jet covered most of 
the estimated values for classical GRBs (from 0.04 rad 
to 0.5 rad), while 9obe varied from 0.0 rad to 1.5 rad 
when an at least 3-fold enhancement of 9jet was exam¬ 
ined. Figure 0] shows the best-fit model functions that 
describe the achromatic rebrightening, i?/i?c-band tem¬ 
poral evolution, and radio brightness. The derived burst 
parameters are also presented in Table [2] There are two 
notable features. The first is the jet opening half angle 
9jet of 0.1 rad, which is cons istent with those of classical 
hard GRBs (see Figure 6 in iFong et al.l 1201^ . The sec¬ 
ond is the large observing angle {9obs "^0.21 rad), which 
corresponds to 9obs ‘29jet- 

To compare with other samples such as an off-axis 
origin XRF candidate, classical hard GRB, XRR, and 
the on-axis orphan GRB candidates, we considered 
080330, 990510, 120326A, 131030A, and PTFllagg (Ta¬ 
ble [2]). The burst parameters for these bursts, except 
080330, were also determined by using the same box 
ht code on the basis of multifrequency afterglow obser¬ 
vations (I van Eer ten et al.ll2012t lUrata et al.ll20l4 120151 : 
ICenko et al.ll2013h . To compare the burst parameters of 
XRF080330 with those of XRF020903 and others, we also 
performed the multiband light curve modeling with the 
sam e code by using availa ble multiband photometric re¬ 
sults [CuidorzreF^ (|20^. XRF080330 was detected by 
Swift and an afterglow in optical and IR bands showed 
the achromatic brightening phase that could be inter - 
preted using the off-axis jet model (|Guidorzi et al.l[200^ . 
However, was not well characterized, as being less 

than 88 keV (IGuidorzi et al.ll^OOOf l. because of the rel¬ 
atively high and narrow energy band of Swift/JiAT. In 
addition, the large uncertainties in the of fluence ratio of 
the prompt emission (l.bj'jQ'g) were still compatible with 
an XRR category base d on the modified XRF /XRR def¬ 
inition for Swift/BAT (iSakamoto et al.l|2008ll . As shown 
in Figure 01 the achromatic brightening light curves 
were well modeled, and the derived burst parameters 
are presented in Table [2] One of the key features was 
^ohs(=0.12 rad), which was equal to 0jet(O.12 rad). Thus, 
the achromatic brightening is explained by the off-axis 
jet model. However, our viewi ng angle (9nhs ^ 9jet) 
was not larger as estimated by IGuidorzi et al.l (|2009ll 
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~ 1.5-2 X djet). 

This enables a comparison of the burst parameters, 
which are presented in Table [2j We also list the observed 
features of the prompt emission. The observed values 
are widely distributed to represent the variety of GRB 
classes. By contrast, the parameters obtained through 
the afterglow modeling are comparable to each other, 
except for Oohs- Thus, the off-axis jet model is suitable 
for explaining the diverse afterglow light curves and the 
GRB category. 


5.2. Small values of A'peafc ^iso with the Off-axis 

model 

Using the values of Oj^t and 0obs obtained by fit¬ 
ting afterglow light curves (Table [U, we discuss on 
observed small values of and Ei^o to verify 

the off-axis jet model. We adopted a simple model 
with a top-h at profile of t he prompt emission of rel¬ 
ativistic let (IGranot et al.l 119991: iWoods fc LoebI 119991 : 
lloka fc Nakamura! 120011 : lYamazak i et al.l 1 2002 ). Fo llow¬ 
ing the formalis m derived bv iGraziani et alH 20061) and 
iDonaghvI (1200611 . the peak energy E^ff^f.{6obs) and the 
isotropic energy Eiso(Qobs) were analytically derived as 
functions of Oobs, and the Lorentz factor of the jet 
7 = (1 — We defined the ratios i?i and R 2 as 


TPSTC (n \ 

Epeak^Pobs) 

~ 

_ 2(1 - /3)(1 - fcosOjet) 

2 - /3(1 -I- cos Ojet) 

f{l3 - cos9obs) - fiflcosdjet - cos9obs) 
g{l3 - cosdobs) - g{/3cos9jet - cos6obs) 

_ Eiso {9obs) 

Eiso{0) 

_ (l-/^)^(l-/^COSdjet)^ 

/3(1 - cos9jet)[2 - /3(1 -I- cos9jet)] 
x[f{/3- cos9obs) - f {Pcos9jet - 9obs)] 


, ( 1 ) 


( 2 ) 


where 


f{z) = 
7^(272 


9iz) = 


l)z^ -I- (372 sin2 9obs - l)z -I- 2 cos 9obs sin^ 9obs 

1^2 _|_^-2 gij^2 

272 z -I- 2 COS 9 obs 
2 -|- 7“2 sin2 9obs\i 


(3) 

(4) 


For fixed 9jet a nd 7 , it can be seen that for 9obs > 9jet 
(see Figure 2 of lDonagh\ll2006ll . both and Eiso de¬ 

creases with an increase in 9obs because of the relativistic 
beaming effect. 

We consider the case of XRF 020903. In the following, 
we fix 9obs = 0.21 rad and 9jot = 0.1 rad. Subsequently, 
ratios Ri and R 2 are calculated for given 7 . For exam¬ 
ple, we find Ri = 3.45 x 10“^ and R 2 = 1.23 x 10“® if 
7 = 100. Glassical hard GRBs typically have = 

500 keV (INava et al.ll2012ll . and therefore we first assume 
= 500 keV. We then find that = 

0.21) = 4.2 keV, which is close to the observed result 
for XRF 020903. Furthermore, for the observed value 


of Eiso, we consider Eiso{9obs = 0.21) = 1.4 x lO'^^erg, 
resulting in Eiso{0) = 1.9 x 10^'^erg. In summary, if a jet 
with 9jet = 0.1 and 7 = 100 is seen on-axis {9obs = 0), we 
would have = 500 keV and Eiso = 1.9 x lO^^erg, 

which is almost consistent with the E^f^^ — Eiso relation 
(|Amati et al. 11200^ INava et "3111201211 . The empirical re¬ 
lation could be an indicator of GRB category although 
the background physics are not yet fully understood. In 
fact, classical hard GRBs show the relation, but short 
GRBs do not exhibit it. Hence, the observed small val¬ 
ues of and Eiso for XRF 020903 are naturally ex¬ 
plained by the off-axis jet model. 

We also discuss parameter dependence. In the fol¬ 
lowing, we fix 9obs = 0.21 rad, 9jet = 0.1 rad, and 
Eiso{9obs) = 1-4 X 10^®erg. The observed best ht value 
Epeaki^obs = 0.21) = 3.3 keV is reproduced for 7 = 112 
and Eplf^f.{0) = 500 keV, and we obtain Eiso{0) = 
3.0 X 10®"^erg. For £’*((((^(0) = 250 keV, we need a smaller 
7 value (79) is required to obtain E^Jf^f,{9obs = 0.21) = 
3.3 keV, and we hnd Uiso(O) = 7.6 x lO^^erg. Similarly, 
when we consider U®g^j,(0) = 1 MeV, we have 7 = 160 

and Uiso(O) = 1.3 x lO^^erg. For these cases, the on- 
axis values, Flpeafc(O) and Eiso{0) remain within the do- 
scatter of the i^peak “ ^iso relation. These values of 7 
are consistent with those measured u sing the onset of 
optical afterglows. iLiang etldl (|201fl[ l reported 7 mea¬ 
surements that were distributed from ^ 90 to ~ 600 
as well as a significant proportion of events that exhib¬ 
ited values below 200. However, the values are lower 
than the predi ction (7 > 300 for all GRBs) made by 
lDonagh\^ (|2006f) for reducing the number of unseen events 
away from the observ ed - Eiso and Epl^j. - E^ 

(jGhirlanda et al.l l2004fl relations. Here, E..^ is the jet 
collimation-correc ted energy. The simplest off-axis jet 
models for XRFs (|Yamazaki et al.l[2002Ll200^ adopt low 
values of 7 (^ 100 ) and predict a large number of events 
away from these rela t ions t hat are not observed (see Fig¬ 
ures 4-9 in IDonaghvI (|2006[B . Therefore something more 
complicat ed procedures m ust be required. For one of 
examples, IDonaghvI (|2006f ) described that a inverse cor¬ 
relation between 7 and the opening solid angle of the 
GRB jet has the effect of greatly reducing the visibility 
of off-axis events.Another poss ibility is to conside r more 
complicated jet structure ('e.g. iSalafia et al.ll2ni5h . 

The spectral peak energy of XRF020903 is conser- 
vatively determ i ned o nly the upper limit, reported by 
iSakamoto et al.l (|2004ll . Although XRF020903 also fol¬ 
lowed the — Eiso relation when we employed the 

yalue estimated by the constrained Band function with 
the 90% confident leyel, four of the fiye HETE-2 XRFs 
(including XRF 020903) with the lowest upper limits on 
Epllk were potential outlier eyents of the Epl'^^ - Eiso 
and the Liso — relations. Because these eyents are 
outliers in the (Figure 1 in ISakamoto et al.l 

()200,5h and Figure 9 in iLamb et'm (l2005h l and the 
^peak ~ Se planes (Figure 2 in IDonaghvI 120061 1. Here, 
Ffi and Se denote the peak photon number flux of the 
burst and the energy fluence, respectively. These results 
indicate that the nature of these four XRFs — and there¬ 
fore the nature of XRFs in general — might differ con- 
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TABLE 2 

Summary of burst parameters obtained by observations and numerical modeling 


Parameters 

020903 

080330 

990510 

131030A 

120326A 

PTFllagg 

Category 

(keV) 

Eiso (erg) 
z 

XRF 

q q + 1.8 
'^•'^-1.0 

i.4tjy X 10^9 
6.251 

XRF(XRR?) 

< 88 

< 2.2 X 10®2 
1.51 

GRB 

423t« 

2.lt°-3 X 10®3 

1.619 

GRB 

406 ± 22 
3 . 0/32 X 10®3 
1.293 

XRR 
107 . 8/(5 3 

3.2/° ( X 1052 
1.798 

on-axis orphan (?) 

0.5< .s <3.0 

Ojet (rad) 

0.10 

0.12 

0.075 

0.15 

0.14 

0.20 

E (erg) 

5.9 X 1059 

2.3 X 10®9 

1.8 X 10®3 

3.4 X 10®2 

3.9 X 1052 

9 X 1052 

n (cm 

1.1 

9.0 

0.03 

0.3 

1.0 

0.001 

Oaba (rad) 

0.21 

0.12 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0.19 

P 

2.8 

2.1 

2.28 

2.1 

2.5 (fixed) 

3.0 


1.4 X 10-3 

1.6 X 10-1 

4.6 X 10-3 

4.4 X 10-2 

1.0 X 10-3 

4 X 10-2 

ee 

2.9 X 10-1 

1.4 X 10-1 

3.7 X 10-1 

2.7 X 10-1 

6.9 X 10-1 

2 X 10-1 

x'^/dof 

90.9/9 (10.1) 

512.5/125 (4.1) 

1267.2/198 (6.4) 

- 

28.0/20 (1.4) 

- 

Data 

Opt, Radio 

Opt 

X,Opt,Radio 

Opt, ALMA 

Opt 

Opt, Radio 

Ref. 

This work 

This work 

van Eerten et al. (2012) 

Urata et al. f2015) 

Urata et al. 12014) 

Cenko et al. 12013) 


siderably from those of the rest of the XRFs, XRRs, and 
GRBs. However, our result on one of the four outlier 
events contradicts this suggestion. Hence, we may be 
able to explain the remaining three lowest events 

in an identical manner, although additional adjustments 
on the jet geometry or the initial Lorentz factor could be 
required. 

5.3. Prospect with optical surveys 

Finding off-axis orphan afterglows through untargeted 
optical surveys is also a crucial method for establish¬ 
ing a unified picture of GRBs. Several optical time- 
domain surveys performed using telescopes with diam¬ 
eters in the range of l-2m (e.g. iPTF, Pan-STARRSl) 
have reported intriguing new discoveries related to stel¬ 
lar explosions. However, the sensitivities of these surveys 
were insufficient for detecting faint off-axis orphan after¬ 
glows. Hence, off-axis orphan GRB afterglows are yet 
to be observationally confirmed. However, a new wide- 

field-of-view camera-Hyper-Suprime-Gam (HSG)- 

attached to the Subaru telescope and the planned LSST 
have considerable capabilities to detect the first off-axis 
orphan GRB afterglow in untargeted time-domain sur¬ 
veys. 

One of the challenges in generic transient surveys is 
determining candidates from the various types of opti¬ 
cal transients, because the occurrence of GRB orphan 
afterglows is rarer compared with that of known types of 
supernova. These candidate selection can be performed 
using optical photometric survey data and a proper pho¬ 
tometric transient classification. The photometric clas¬ 
sification involves seven steps: (1) finding transient can¬ 
didates by generating differential images, (2) generat¬ 
ing light curves for transient components, (3) identify¬ 
ing host galaxies, (4) determining a transient location 
in their hosts, (5) matching the candidates with known 
sources in various catalogs, (6) matching light curves and 
color evolution, and (7) estimating the photometric red- 
shift of ho sts. For Step (F) and ( 2). we can employ the 
algorithm (lAlard fc Luoton 1119981) as described in §4 and 
lUrata et al.l (I2012bl) . For'^ep (3), (4), and (7), a con¬ 
siderable number of GRB host galaxy observations have 
been performed in the optical and near-infrared range. 
Systematic unbiase d observations hav e also been per¬ 
formed using VLT (jHiorth et al.ll20l^ . The brightness 
range of host galaxies for GRBs for z < 2 is 23.0-26.5 


mag fe.g.. lBergei<2010ll . Since the redshift range of a con¬ 
siderable number of orphan GRB afterglows with HSG 
surveys extends up to z ^ 2, HSG images (e.g. ref¬ 
erence images for PSF matched subtraction) are suffi¬ 
ciently deep to detect these host galaxies. Thus, we can 
also perform the photometric redshift for the host galax¬ 
ies. Photome tric redshift for GRB ho st galaxies is also 
effective fe.g.. IGhristensen et al.ll20bll . Hence, the light 
curve expectations basis of off-axis jet origin of XRFs are 
crucial to establish a effective candidate searches. 

Figure [7] shows the expected off-axis afterglow light 
curves at z = 1 along with various observing angles 
(Qobs = Ojet, ‘20jet, and 39jet) in the optical r-band, 
ALMA Band3, and JVLA 8.5 GHz, obtained from the 
light curve modeling of XRF020903. The orange dashed 
lines in Figure [7] indicate the sensitivity limit for 1 h 
exposure for each instrument. Hence, off-axis orphan 
GRB afterglows (up to ^ 39jet) can be detected in op¬ 
tical time-domain surveys by using 8m class telescopes. 
Follow-up radio observations are also crucial for confir¬ 
mation of orphan afterglows and identification of con¬ 
straints on their physical parameters, as shown in Fig¬ 
ure [71 Because radio temporal evolution is substantially 
slower than that of optical temporal evolution, long-term 
monitoring by using ALMA and/or JVLA with reason¬ 
able exposure (^Ih) requires the confirmation of the op¬ 
tical candidates. 

6. CONCLUSION 

We studied XRFs on the basis of redshift measure¬ 
ments, multifrequency afterglow modeling, and spectral 
peak measurements (A®™^) to verify the off-axis jet 
model and to provide feedback to ongoing and planned 
optical time-domain surveys, which have considerable 
potential for detecting off-axis orphan GRB afterglows. 
Because off-axis orphan GRB afterglows are produced as 
a natural consequence of GRB jets production, the con¬ 
firmation of the off-axis origin of XRFs is the necessary 
to conduct off-axis orphan GRB afterglow surveys. 

To verify the off-axis jet model, we selected XRF020903 
by considering the three aforementioned sample selection 
factors. For this event, we reduced the archived data of 
Subaru to describe the optical light curves, and found 
achromatic rebrightening at ^0.7 days. Using these op¬ 
tical results and radio data obtained from literature, we 
performed afterglow light curve modeling with the boxfit 
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Fig. 7. — Expected off-axis GRB afterglow light curves with observing angle of 26jf,t (red solid line) and 36jet (blue solid line) at ^ = 1 
in the r'-band (top), ALMA Bands (middle), and JVLA 8.5 GHz (bottom) obtained from the off-axis jet modeling of XRF020903. The 
dash orange lines indicate the sensitivity limit of Subaru (r'-band), ALMA (Bands), or JVLA (8.5GHz) with f h exposure. 


code and found that the off-axis jet model {Oobs ^ ‘^Qjet) 
could explain the achromatic rebrightening, R/Rc-hawA 
temporal evolution, and radio brightness. 

We also compared the burst parameters of XRF020903 
with those of other categories of events, such as a classical 
hard GRB, an XRR, and an on-axis orphan GRB candi¬ 
date. For XRF080330, we performed light curve model¬ 
ing in a manner similar to that used for XRF020903 by 
using optical data from literature, and we confirmed that 
the off-axis jet model (0obs ~ Ojet) could describe the op¬ 


tical afterglow light curves. We also listed the observed 
features of the prompt emission for each event. The ob¬ 
served values were too widely distributed to represent the 
classical hard GRB, XRR, and XRF. By contrast, the 
parameters obtained from the afterglow modeling were 
comparable to each other, except for Oobs- Thus, the off- 
axis jet model was found to be suitable for explaining the 
diverse afterglow light curves and the GRB category. 

We also verified the observed small values of and 
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Eiso by adopting a simple model with a top-hat profile 
of the prompt emission of a relativistic jet. The parame¬ 
ters and Eiso{0obs) were analytically derived 

as fnnctions of 9obs, 9jet, and 7 = (1-/3^)“^/^. By fixing 
Oobs and Ojet as 0.21 and 0.1 rad, respectively, we evalu¬ 
ated 7 , E^l2k^ Eiso observed from the on-axis of the 
jet {Epl^i.{0) and Eiso{0))- These expected values were 
consistent with those of classical hard GRBs, and the 
observed small values of and Eiso of XRF 020903 
could be naturally explained by the off-axis jet model. 

Finally, we expected off-axis orphan GRB afterglow 
light curves at 2 ; = 1 along with three viewing angles 
on the basis of the XRF afterglow light curve model¬ 
ing. To detect these light curves, especially afterglows 
with a larger viewing angle {9obs > an 8 -m class 

telescope with wide-held imagers, such as the LSST and 
Subaru/HSC, is required. Off-axis orphan GRB after¬ 
glows up to ^ 30jet can be discovered by performing 
time-domain surveys with an 8 -m class telescope. Be¬ 
cause such optical time-domain surveys also detect nu¬ 
merous other optical transients, we presented expected 
radio afterglow light curves for the conhrmation and de¬ 
termination of burst parameters. Radio light curves can 
be monitored using ALMA and JVLA with reasonable 


Alard & Lupton 1998, ApJ, 503, 325 
Amati, L., 2002, A&A, 390, 81 

Baba, H., et al., 2002, ASP Conference Series, Vol. 281, 298 
Band, D., 1993, ApJ, 413, 281 

Barthelmy, S. D., Barbier, L. M., Cummings, J. R., et al. 2005, 
Space Sci. Rev., 120, 143 
Berger, E. 2010, ApJ, 722, 1946 

Bersier, D., Fruchter, A. S., Strolger, L.-G., et al. 2006, ApJ, 643, 
284 

Bertin, E., Mellier, Y., Radovich, M., et al. 2002, Astronomical 
Data Analysis Software and Systems XI, 281, 228 
Cenko, S. B., Kulkarni, S. R., Horesh, A., et al. 2013, ApJ, 769, 
130 

Christensen, L., Hjorth, J., & Gorosabel, J. 2004, A&A, 425, 913 
Donaghy, T. Q. 2006, ApJ, 645, 436 

Fong, W., Berger, E., Margutti, R., et al. 2012, ApJ, 756, 189 
Ghirlanda, G., Ghisellini, G., & Lazzati, D. 2004a, ApJ, 616, 331 
Granot, J., Piran, T., & Sari, R. 1999, ApJ, 513, 679 
Graziani, C., Lamb, D. Q., & Donaghy, T. 2006, Gamma-Ray 
Bursts in the Swift Era, 836, 117 (astro-ph/0505623) 

Guidorzi, C., Clemens, C., Kobayashi, S., et al. 2009, A&A, 499, 
439 

Heise, J., Zand, J. I., Kippen, R. M., & Woods, P. M. 2001, 
Gamma-ray Bursts in the Afterglow Era, 16 
Heise, J. 2003, Gamma-Ray Burst and Afterglow Astronomy 
2001: A Workshop Gelebrating the First Year of the HETE 
Mission, 662, 229 

Henden, A. 2002, GRB Goordinates Network, 1571, 1 
Hjorth, J., Malesani, D., Jakobsson, P., et al. 2012, ApJ, 756, 187 
Huang, Y. F., Dai, Z. G., & Lu, T. 2002, MNRAS, 332, 735 
loka, K., & Nakamura, T. 2001, ApJ, 554, L163 
Jager, R., Mels, W. A., Brinkman, A. G., et al. 1997, A&AS, 125, 
557 

Kaneko, Y., Preece, R. D., Briggs, M. S., et al. 2006, ApJS, 166, 
298 

Kippen, R. M., Woods, P. M., Heise, J., et al. 2001, Gamma-ray 
Bursts in the Afterglow Era, 22 

Lamb, D. Q., Donaghy, T. Q., & Graziani, C. 2005, ApJ, 620, 355 
Liang, E.-W., Yi, S.-X., Zhang, J., et al. 2010, ApJ, 725, 2209 
Lien, A., et al. ApJS in prep. 

Magnier, E. A., Schlafly, E., Finkbeiner, D., et al. 2013, ApJS, 
205, 20 

Magnier, E. 2006, The Advanced Maui Optical and Space 
Surveillance Technologies Gonference, 


exposure. 


This work is partly supported by the Ministry of 
Science and Technology of Taiwan grants NSG 100- 
2112-M-008-007-MY3 and MOST 103-2112-M-008-021- 
(YU). Subaru Suprime-Gam data were acquired through 
SMOKA, which is operated by the Astronomy Data Cen¬ 
ter, National Astronomical Observatory of Japan. The 
PSl Surveys have been made possible through contri¬ 
butions of the Institute for Astronomy, the University 
of Hawaii, the Pan-STARRS Project Office, the Max- 
Planck Society and its participating institutes, the Max 
Planck Institute for Astronomy, Heidelberg and the Max 
Planck Institute for Extraterrestrial Physics, Garching, 
The Johns Hopkins University, Durham University, the 
University of Edinburgh, Queen’s University Belfast, the 
Harvard-Smithsonian Center for Astrophysics, and the 
Las Cumbres Observatory Global Telescope Network, In¬ 
corporated, the National Central University of Taiwan, 
and the National Aeronautics and Space Administra¬ 
tion under Grant No. NNX08AR22G issued through the 
Planetary Science Division of the NASA Science Mission 
Directorate. 


Malesani, D., Fynbo, J. P. U., Jakobsson, P., Vreeswijk, P. M., Sz 
Niemi, S.-M. 2008, GRB Coordinates Network, 7544, 1 
Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999 
Miyazaki, S., et al. 2002, PASJ, 54, 833 

Miyazaki, S., Komiyama, Y., Nakaya, H., et al. 2012, Proc. SPIE, 
8446, 84460Z 

Nava, L. et al. 2012, MNRAS, 421, 1256 

Ohno, M., Fukazawa, Y., Takahashi, T., et al. 2008, PASJ, 60, 361 
Ouchi, M., Shimasaku, K., Okamura, S., et al. 2004, ApJ, 611, 660 
Pierfederici, F., Valdes, F., Smith, C., Hiriart, R., &; Miller, M. 
2004, Astronomical Data Analysis Software and Systems 
(ADASS) XIII, 314, 476 

Pierfederici, F. 2006, Astronomical Data Analysis Software and 
Systems XV, 351, 433 
Rhoads, J. E., 1999, ApJ, 525, 737 

Sakamoto, T., Lamb, D. Q., Graziani, C., et al. 2004, ApJ, 602, 
875 

Sakamoto, T., Lamb, D. Q., Kawai, N., et al. 2005, ApJ, 629, 311 
Sakamoto, T., Hullinger, D., Sato, G., et al. 2008, ApJ, 679, 570 
Salafia, O. S., Ghisellini, G., Pescalli, A., Ghirlanda, G., Sz 
Nappo, F. 2015, arXiv:1502.06608 
Sari, R., Piran, T., Sz Halpern, J. P. 1999, ApJ, 519, L17 
Schlafly, E. F., Finkbeiner, D. P., Juric, M., et al. 2012, ApJ, 756, 
158 

Schlafly, E. F., &; Finkbeiner, D. P. 2011, ApJ, 737, 103 
Serino, M., Sakamoto, T., Kawai, N., et al. 2014, PASJ, 66, 87 
Shirasaki, Y., Kawai, N., Yoshida, A., et al. 2003, PASJ, 55, 1033 
Soderberg, A. M., Kulkarni, S. R., Berger, E., et al. 2004, ApJ, 
606, 994 

Tonry, J. L., Stubbs, C. W., Lykke, K. R., et al. 2012, ApJ, 750, 
99 

Urata, Y., Huang, K., Yamaoka, K., Tsai, P. P., Sz Tashiro, M. S. 
2012, ApJ, 748, L4 

Urata, Y., Tsai, P. P., Huang, K., et al. 2012, ApJ, 760, LLll 
Urata, Y., Huang, K., Takahashi, S., et al. 2014, ApJ, 789, 146 
Urata, Y., Huang, K, et al. 2015, ApJ, in prep, 
van Eerten, H., van der Horst, A., & MacFadyen, A. 2012, ApJ, 
749, 44 

von Kienlin, A., Meegan, C. A., Paciesas, W. S., et al. 2014, 
ApJS, 211, 13 

Woods, E., & Loeb, A. 1999, ApJ, 523, 187 

Yamaoka, K.., et al. 2009, PASJ, 61, S35 

Yamazaki, R., loka, K., Sz Nakamura, T. 2002, ApJ, 571, L31 

Yamazaki, R., loka, K., and Nakamura, T., 2003, ApJ, 593, 941 


Off-Axis GRB and Orphan GRB afterglow 


11 


Zhang, B., Dai, X., Lloyd-Ronning, N. M., & Meszaros, P. 2004, 
ApJ, 601, L119 



